CPT-P09-2006 



UNE INTERSECTION DE QUADRIQUES LIEE A LA SUITE DE STURM 

Oleg Ogievetsky^, Vadim Schechtman^ 



TABLE DES MATIERES 

Première Partie. Formules 

§ 1. Introduction 2 

§ 2. Algèbre S 4 

§ 3. Début de la démonstration du théorème 1.5 8 

§ 4. Formule (^4) 10 

§ 5. Formule (S) 11 

Deuxième Partie. Polynômes d'Euler et déterminant de Cauchy 

§ 1. Nombres fi{j)i 14 

§ 2. Polynômes d'Euler et fonction hypergéométrique 16 

§ 3. Asymptotiques 17 

Bibliographie 20 



'^Centre de Physique Théorique, Luminy, 13288 Marseille, France (Unité Mixte de Recherche 6207 du CNRS et des 
Universités Aix-Marseille I, Aix-Marseille II et du Sud Toulon - Var; laboratoire affilié à la FRUMAM, FR 2291) et 
Institut de Physique P.N. Lebedev, Département Théorique, Leninsky prospekt 53, 119991 Moscou, Russie; 
email: oleg@cpt.univ-mrs.fr 

^Laboratoire Emile Picard, UFR MIG, Université Paul Sabatier, 31062 Toulouse, France; 
email: schechtman@math.ups-tlse.fr 



1 



PREMIÈRE PARTIE. FORMULES 



§ 1. Introduction 

1.1. Cet article est une variation sur un thème de [Jacobi]. 
Soit 

f{x) = a„a;" + an-ix"'~^ + . . . + oq 

un polynôme de degré n > à coefficients dans un corps de base É de caractéristique 0. Rappelons que 
la suite de Sturm de /, 

f = ifo, fl,f2, ■ ■ ■) , 

est définie par récurrence : on pose fo{x) = f{x), fi{x) = f'{x) et pour j > 1 fj+i est le reste de la 
division euclidienne de par fj, avec le signe opposé : 



fj-i{x) = qj-i{x)fj{x) - fj+i{x), 



(1.1.1) 



deg/j+i(x) < degfj{x), cf. le célèbre mémoire [Sturm]. 

Dans cette note on propose des formules explicites pour les coefficients des polynômes fj en termes 
des coefficients de /. Plus gcncralcmcnt, on donnera des formules analogues pour les membres de 
l'algorithme d'Euclidc correspondant à deux polynômes quelconques /i,/2 de degrés n — l,n — 2. 

Notre point de départ est une algèbre 05, quotient de l'anneau de polynômes en variables b{i)j {i > 
1) j > 2z) par certains rélations quadratiques, cf. (1.7.1) ci-dessous. Nos formules sont des conséquences 
des identités dans *B, analogues des rélations de Pliicker. 

1.2. Pour énoncer le résultat, introduisons les quantités quadratiques 

i-i 

p=0 

i > 1, î > 2_7- Ici on pose Oj = pour z < 0. Par exemple, 

6(l)î = nianan-i - {n - i + l)a„_ia„_î+i . 



1.3. Ensuite on introduit, pour m > 2, les matrices (m — 1) x (m — 1) symétriques 



C{m) 



De plus, pour i > on définit une matrice "décalée" C{m)i : elle est obtenue en remplaçant dans 
C{m) la dernière ligne par 

( &(l)m+i 6(2)m+i+l b{^)m+i+2 b{4:)m+i+3 ■ ■ ■ ^'(m - l)2m+i-2 ) • 



/ Kl)2 


6(1)3 


6(1)4 


6(1)5 . 


6(1)™ 


\ 


6(1)3 


6(2)4 


6(2)5 


6(2)6 . 


6(2)^+1 




6(1)4 


6(2)5 


6(3)6 


6(3)7 . 


6(3)^+2 




6(1)5 


6(2)6 


6(3)7 


6(4)8 . 


6(4)^+3 




V b{l)m 


6(2)^+1 


6(3)m+2 


6(4),n+3 • 


. 6(m - l)2m-2 


/ 
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Donc C(m)o = C{m). On pose 

c{m)i := det C(m)i, c(m) := c(m)o 

En particulier, 

C(2), = h{l)i+2 

Il est commode de poser 

{n-i)an-i 

c(l)i := 



i > 0, c(l) := c(l)o = 1. 

1.4. Puis on définit les nombres 7j, j >1 par récurrence : 

1 c(j - 1)' 

j > 2. Autrement dit, 

i=l 

oii = nan si j est impair et l/(n^a„) sinon. 

Les nombres 71, . . . , 7^ sont donc bien définis si tous les nombres c(2), c(3), . . . ,c{j—l) sont différents 
de zéro. 

1.5. Théorème. Supposons que deg/j = n — j, donc deg/j = n — i pour z < j. 
Alors pour tous z < on a c{i) 7^ et 

p=0 

En particulier, le coefficient dominant de fi{x) est égal à jic{i). 

1.6. On vérifie aussitôt que 

b{k)i - b{k - l)i = c(l)fc_i6(l),_fc+i - c(l),_fc&(l)fc (1-6.1) 
pour tous k >2, i>2k — 2. Par exemple, 

b{2)i - b{l)i = c(l)i6(lX_i - c(l)i_26(l)2, 

6(3), - b{2)i = C(l)26(l)i_2 - C(l)i_36(l)3, 

etc. Il s'en suit que tous les b(j)i, j > 2, sont expressibles en termes de c(l)p et c(2)p = 6(l)p_|_2, p > 0. 

1.7. Les formules (1.6.1) impliquent que les nombres b{i)j satisfont aux relations quadratiques 
suivantes : 

ibik)i - b{k - ■ 6(1),- 

(1.7.1) 

= {b{j)i-k+j - b{j - l)î-fe+j) • ^(l)fc - {Kj)k+j-l - Kj - l)fc+j-l) • &(l)i-fe+l 
On verra que la preuve de 1.5 ne dépend que des relations (1.7.1). 



On formalise la situation en introduisant une algèbre quadratique correspondante, cf. § 2 ci-dessous. 

1.8. Maintenant soient 

= aox"-'^ + aix"~^ + . . . 

et 

/2(x)=/3ox"-2 + /3ix"-=^ + ... 

deux polynômes arbitraires de degrés n — 1, n — 2. On définit /j, j > 3 à partir de /i, J2 par les formules 
de l'algorithme d'Euclide (1.1.1). 

Posons 

c{l)i :=^, b{l)i+2 -Pi, i>0. 
«0 

Définissons les nombres b{k)i, k >2 par récurrence sur k, à partir des formules (1.6.1). 
Définissons les nombres c{m)i, m > 2, par les formules 1.3. 
Enfin, on pose : 

71 = "0 , 72 = 1 , Ij+i = 7j-i— ^^^p— 

Alors on a 

1.9. Théorème. Supposons que deg/j = n — j, d'où, deg/j = n — i pour i < j. 
Alors pour tous i < j, on a c{i) 7^ et 

n—i 
p=Q 

En particulier, le coefficient dominant de fi{x) est égal à 7ic(i). 
Cf. [Jacobi], section 15. 

1.10. Dans la Deuxième Partie on présente un exemple numérique. Là, les déterminants de Cauchy 
apparaissent dans les asymptotiques des coefficients dominants de la suite de Sturm pour les polynômes 
d'Euler. 



§ 2. Algèbre <B 

2.1. On peut réécrire les relations (1.7.1) sous la forme suivante : 

det( ]-det( ^^^^^ 6(j-l),+.-i\ 

''''' [ b{j - l),+,_fe b{k - 1), ) 1^ 6(l),_fe+i b{k)i J 

+ d^*f tm^'"' ] =HK3)i-^'{k,j),+A"{k,j), = Q 

V b(l)i-k+i b[j)i+j_k I 



(2.1.1) 



2.2. On définit une algèbre quadratique *8 comme une ï-algèbre commutative engendrée par les 
lettres b{i)j, i,j G Z, modulo les relations (2.1.1), oii i,j,k G Z. 

(D'ailleurs, dans tout le paragraphe qui suit on peut remplacer le corps de base î par un anneau 
commutatif quelconque.) 
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2.3. Le but de ce paragraphe est d'écrire certaines relations entre les déterminants n x n dans 05 
qui généralisent (2.1.1). 

On fixe un nombre entier n>2. Soient mi, . . . , m„, i des entiers. 

On définit 2n + 2 vecteurs Vj,Wj E f^, j = 1, . . . ,n + 1 : 

— m„+l jj 

1 < J < (suivant l'usage, x signifie que l'on omet la composante x). 



Puis 



Vi = {h{mi - l)j+m,i_m,„,6(m2 - l)i+m2_m„,...,6(m„_i - l)i+„„_-^_^„,6(mn - l)i) , 

V2 = {b{mi - \)rm^mr,-U K'fn2 - l)m2+m„-l, ■ ■ • , - l)m„_i+mn-l, ! 

^3 = (6(mi - l)mi+mn-i-l, b{^2 - l)m2+m„_i-l, • • • , &(m„_2 - l)m„_2+m„_i-l, 

^('Tîn— l)m„-i+m„— 1) ^('Tln— l)î+m„-i— m„) j 
t;4 = (6(mi - l)mi+m„-2-l^ h{m2 - l)m2+m„_2-l) • • • ) b{mn-2, " l)m„_3+m„_2-l) 
6(mn_2)m„-2+m„-i-l, ?'(ï^^n-2)m„-2+m„-l, &('Tln-2)î+m„-2-m„) , 



= (6(mi - 1) &("l2)m2+m3-l, fc("^2)m2+m4-l, • • • , K"^2)m2+m„-l, b{m2)i 

+1712— mn J 

Vn+1 = (&("ll)mi+m2-l) ^("T'l)mi+m3-l! • • • ) ^("T'l)mi+mn-l) &("ll)i+mi-m„) • 



2.4. Soit 



Xll 



M 



Xn-2,1 ■ ■ ■ Xn-2,n+l 



une matrice (n — 2) x (n + 1) sur 03 ; soit Mi, i = 1, . . . , n + 1, ses sous-matrices (n — 2) x n. Pour 
écrire Mj, on enlève donc la î-ième colonne de M. 

Maintenant on va définir n + 1 matrices n x n 

Dj = Dj{mi, . . . ,m„;M„+2-j)i, 

j = 1, . . . ,n + 1. On pose : 



Mn+l 
\ VI J 



j = 2, . . . ,n + 1. Ici (.)* désigne la matrice transposée. 
Enfin, on pose 



Aj = Aj{mi, . . . ,mn; Mn+2-j)i = det Dj{mi, ... ,mn; Mn+2- 



j = l,...,n + l. 



5 



Considérons la somme alternée 



n+l 



R{n;mi,...,mn;M)i = ^ (-iy+^Aj(mi, . . . ,mn;M„+2-j)i • 

2.5. Exemple, n = 2. Dans ce cas il n'y a pas de matrice M ; trois nombres entiers sont donnés : 
mi,m2 et i. On aura 6 vecteurs : 

m = {b{l)mi,b{l)m2), W2 = {h{l)mi,h{l)i-m2+\), m = (&(l)m2 , ^(l)i-m2+l) 



et 



Il s'ensuit : 



Vi = {b{mi - l)i+rni-m2,K^2 - V2 = {b{mi - l)rai+m2-l, K^2)i), 

V2 = {b{mi)mi+m2-l,b{'^l)i+mi~m2) ■ 



R(2-mi 7712)- = det ( '^^^^'^^ ^^^^"^ 



b{l)mi b{mi - 1) 



b{l)m2 &("î-l)mi+m2-l 



\^ 6(l)i_m2+l b{m2)i j y b{m2)i+mi-m2 ) 

On reconnaît là la partie gauche de (2.1.1) pour (j, k) = (mi, 777,2). H en découle que R{2\ mi, 1712)1 = 0. 
2.6. Exemple, n = 3. Dans ce cas la matrice M se réduit à 4 éléments : 

M = Xi X2 X3 X4^ . 
L'expression R{3; 777i, 7772, 7773; M)j prend la forme 



R{S; 7771, 7772, 1^3', M)j = det 



Xi 



b(.l)m2 
X2 



X3 



\ 6(7771 - l)î+mi-m3 b(m2 - I)i+m2-m3 ^("^3 - l)i / 



/ 6(l)mi Xl 6(7771 - l)^^+m3_i 



-det 



6(1)^1 Xl 6(7771 - l)mi+m2-l \ 



6(l)m2 X 2 6 ( 7772 " I)m2+m3-l + det 6(1)^3 X3 6(7772)fn2+m3- 



6(l)i_m,3+l X4 6(7773)4 



— ms+l 3^4 6(7772)i+m2— ma / 



6(1) 

m2 -^2 6(777i)mi+m2-l 
-det| 6(1)^3 X3 6(777i)mi+m3-l 

b{l)i —ma 

Calculons cette expression. 

On développe le premier déterminant suivant la deuxième ligne et les autres suivant les deuxièmes 
colonnes : 

Ai(3; 7771, 7772, 7773; M4)j = -Xl Al (2; 7772 , 7773)1 + X2 Al (2; 777i , 7773)1 - a;3Ai(2; 777i, 7772)i+m2-m3 , 
A2(3; 7771, 7772, 7773; M3)j = -XiA2(2; 7772, 7773)i + X2A2(2; 777i, 7773), - X4A1 (2; 777i , 7773)^2+m3-l • 

Puis 



A3(3; 7771, 7772, ^773; M2)i = -XiA3(2; 7772 , 7773)i + a;3A2(2; 7771, 7772) i+m2-m3 " a;4A2(2; 7771, 7773)r„2+rn3-l 
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et 



A4(3;mi,m2,m3;Mi)j = -X2A3(2; mi, 771,3)1 + X3A3(2; mi, m2)î+m2-m3 - X4.As{2;mi,m3)m2+m3-i 
Pour abréger les notations on introduit des vecteurs entiers : 

(il, i2, is, ù) := {i, i,i + m2- 1713, 777,2 + 777.3-1) , 

fl = (7771, 7772, ma) , 
fJ-l = (^,2,7773), fj,2 = (7771,7773), /X3 = (7771,7772) . 

On peut réécire les formules ci-desssus sous une forme matricielle : 



{ Ai(3;m;M4), \ 
-A2(3;/7;M3)i 
A3(3;/x;M2)i 
V -A4(3;/x;Mi)i / 



/ -Ai(2;^i)i, Ai(2;;U2)^2 -A,{2; fis)^, \ 

A2(2;/7i)i, -A2(2;/72)i2 Ai(2;//2)i4 



-A3(2;/7i) 



A2(2;/X3)i3 -A2(2;/X2)i4 



V 



/^Xi \ 

a;2 



A3(2; 

/^2)î2 ~A3(2;/X3)i3 A3(2;/X2)î4 / 
En rajoutant : 

i?(3; 7771, 7772, 7773; M)i = -xi • | Ai(2; /xi)ii - ^2(2; /xi)ii + ^3(2; Ati)ii | 

+X2 ■ I Ai(2; /72)i2 - ^2(2; ^l2)^2 + A3(2; ^l2)i}^ - X3 ■ |Ai(2; fi^hs - A2(2; ^3)^3 + ^3(2; ^3)^3! 

+a;4- |Ai(2;/72)i4 - A2(2;/72)i4 +A3(2;/72)i4| = . 

Le théorème ci-dessous généralise ces exemples. 
2.7. Théorème. On a 

R{n;mi, . . . ,mn;M)i = 

pour tous 77, 7771, • • • , M et 7. 

Démonstration : elle se fait par récurrence sur 77. Le cas 77 = 2 est l'exemple 2.5. 
Le passage de 77 — 1 à 77 suit l'exemple 2.6. 
Posons pour abréger 

/x = (7771,..., 7n„) . 
À partir de cela, on introduit 77 + 1 vecteurs fXj G Z"~^ : 



j = 1,...,77, et 



Hj := (7771, . . . , 777^, . . . ,777„), 
:= ("71, . . . ,777„_l, 777„) = /7„_l . 



On définit le vecteur 

(71,72, . . ■,in+l) ■= {hh ■ ■ ■ ,7,7 + 777„_l - 777„, 777^-1 + 777^ - 1) G 2"+^ . 



n— 1 fois 
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En développant les déterminants Aj(n; /x, M„_|_2_j)i, 2 < j < n + 1 suivant la deuxième colonne et le 
déterminant Ai(n; /x, M„+i)j suivant la deuxième ligne, on obtient : 

n+l 

R{n;n;M)i = ^ {-iyxjR{n-l;Hj;Mij)i. . 

Ici Mij est la matrice obtenue en enlevant la première ligne et la j-ième colonne de la matrice M. 
Notre assertion en découle immédiatement par récurrence sur n. 

2.8. On aura besoin d'un cas particulier de ces relations. Prenons 

II = {mi,m2, run) = (2, 3, . . . , n + 1) . 

Pour la matrice M, prenons 

/ 6(1)3 6(2)4 6(2)5 ... 6(2)„+2 6(2),+„+i \ 

6(1)4 6(2)5 6(3)5 ... 6(3)„+3 6(3)i+„+2 

\ 6(l)n 6(2)^+1 6(3)„+2 ... 6(n - l)2n-l - l)i+2n-2 / 

Alors le premier déterminant 

Ai(n; /x; Mn+i)i+2n = c{n + 1)^ . 

On pose par définition : 

c(n + !)• := A2(n; ^u; Afn)i+2n , 

c(n + 1) := A3(n; ^; Mn-i)i+2n ■ 

Par contre, si j > 4 on voit que dans le déterminant Aj(n; /x; Mn+2-j)i+2n la dernière colonne est égale 
à la (n — j ' + 3)-ième colonne, d'où 

Ai{n; n; Mn-2)i+2n = ^5{n; ji; Mn-3)i+2n = ... = A„+i(n; ;u; Mi)i+2n = . 
Donc 2.7 entraine 

2.9. Corollaire. Pour tous n > 3 

c{n)i - c{n)[ + c{nyi = . 
§ 3. Début de la démonstration du théorème 1.5 

3.1. On a 

f{x) = Qnx" + a„_ix"~-^ + . . . + ao . 

La dérivée : 

hix) = f'{x) = nanx^-^ + (n - l)a„_ix"-2 + . . . + = na„|x-i + + . . . + ^| 

= 7i(c(l)ox"-i + c(l)ix"-2 + . . . + c(l)„_i) . 
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3.2. Le quotient de la division euclidienne de deux polynômes f{x) et g{x) = a'j^_^x^ ^+a'^_2X^ ^+ 
... est égal à 



X + 



a„_ia„-i-a„_2a„ 



On fait la division euclidienne : 



/ - {x/n + an-l/n''an)f = '""""""'L^""'^""-^ x"-^ + 3na„a„-3-(n-2)a„-ia.-2 ^n-3 + _ _ _ 

= • (6(l)2x"-2 + h{l)^x^-^ + ... + h{l)n) = -72 • (c(2)ox"-2 + c(2)ix"-3 + . . . c(2)„_2) . 

Donc 

n-2 

/2(x)=72-E c(2),a;"-2-V 

î=0 

Cela démontre l'assertion 1.5 pour j = 1, 2, et l'on procède par récurrence par j. 
3.3. On suppose que l'on a déjà trouvé : 

= 7,_i ■ [c(j - + c(i - + . . . + c(j - l)iX^-^+'-' + ...] 

et 

/i(^) = Ij ■ + + ■ ■ ■ + c{j)ix^-^-' + ...]. 

On fait la division euclidienne : 



/.-i(^) - (- 



7jc(i) 



X + 7j-i 



c(j-l)i-%-l^.c(j)i 



1 



7,-i{c(j - 1). - ■ c(i), - 



1=2 
n-j+1 



c{j) 



c{j - 1)1 - ■ c{j). 



c{j)i-i 



c{j) 



■ X 



n-j+l-i _ 



^ ^ |c(i-l),c(jO'-c(i-l)c(i)c(jOi-c(j-l)ic(jOi-i^ 



,=2 

On pose : 

Q{j)i ■■= c{j - l)ic{jf - c{j - l)c{j)c{j)i - c{j - l)icO-)i_ic(j) + c{j - l)c(j)ic(i)i-i (3.3.1) 
Alors on a : 



1=2 



n-j+1 



Il faut montrer que 



ou 



n—j—l n—j+1 

/,-+i(x)=7,+i E c(j + l),a;"-^-^- = 7,-+i E c(j + l),_2x"-^+^- 

i=0 i=0 



c(j-l)^ 



= 7.-1 ■ ^ 
Donc notre théorème est équivalent à l'identité suivante : 

Q{j)i = -cUfcij + i)i . 



(3.3.2) 
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§ 4. Formule {A) 



4.1. Revenons à notre algèbre 03. 
On considère la matrice n x n 



C{n + l)i-2 



( Kl)2 
Kl)3 



6(1)3 
6(2)4 



h{l)n 6(2)n+l 



h{l)n 



6(l)n+l \ 
6(2)n+2 



b{n - l)2n-2 b{n - l)2n-l 
6(n - l)2n+î-3 6(n)2n+î-2 / 



V 6(l)n+i-l 6(2)n+j 

Donc c(n + l)i-2 = det C(n + l)i-2- 

Si l'on désigne par C(n + l)j_2;p,q la matrice C(n + 1)î-2 avec la p-ième ligne et la g-ième colonne 



enlevée, on aura : 



En plus, on a : 



c(n) = det C(n + l)i-2;n,n , 

C(n)i-i = det C(n + l)i_2:n-l,n 



c(n)^' = detC(n + l)._2.„^2,n 

oii c{n)1 a été introduit dans 2.9. 

4.2. Théorème. Pour tous n,i eZ, n > 3, on a la relation suivante dans 03 

c(n — l)ic(n)^ — c(n — l)c(n)c(n)î — c(n — l)ic(n)i_ic(n) + c(n — l)c(n)ic(n)i_i 

= -c(n - l)2c(n + l)i_2 (F) 

On a vu que notre théorème principal 1.5 est une conséquence de {F) : en effet (F) coïncide avec 
la formule (3.3.2) (avec j remplacé par n). 

À son tour, (F) est une conséquence immédiate de deux formules : 



ou bien 



et 



c(n — l)îc(n) — c(n — l)ic(n)î_i = — c(n — l)c{n)" 

c{n — l)ic(n) — c(n — l)ic(n)j_i + c(n — \)c{n)'l = 
{c(n)i + c(n) • c(n) - c(n)ic(n)i_i = c(n - l)c(n + l)i_2 . 



{A) 
{A') 
{B) 



La démonstration de (B) utilise les relations quadratiques entre les lettres h{i)j. Par contre, (A) 
est "élémentaire", en ce sens que cette identité n'utilise pas de relations entre les lettres b{ï)j. 

Pour démontrer (A), on applique le lemme suivant (une variante des relations de Pliicker) : 

4.3. Lemme. {An) Considérons n vecteurs de dimension n — \,Wi = {wn, . . . , Wi^n-i)^ i = 1, ■ ■ ■ ,n. 
A partir d'eux, on définit n vecteurs de dimension n — 2 : Vi = {wn, . . . ,Wi^n-2)- On pose : 

Wi = det(u;i, ...,Wi,..., Wnf , 

Vij := det{vi,...,Vi,...,Vj,...,VnY • 

Alors 

Vn-2,n-l ■ Wn - K-2,n " + K-l,n • = . 
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(Bn) Considérons n vecteurs de dimension n — 2, Vi = {vn, . . . ,Wî,n-2); i = 1, ■ ■ ■ ,n. Considérons 
les mineurs 

Vij := det{vi, û^, , v^f ■ 

Alors pour chaque i < n — 2, 

En effet, en développant Wi par rapport à la dernière colonne, on obtient : (An). 

Par contre, pour vérifier (-B„), considérons la matrice (n — 1) x (n — 2), = Vi. Alors on 
aura Vij = WJ", j = n,n — l,n — 2. D'un autre côté, en développant les mineurs dans {Bn) : 
Vpq, n — 2<p<q<n par rapport à la i-ième ligne, on obtient les mineurs V^, oii y~ est obtenue de 
W en enlevant la dernière colonne. On vérifie que se réduit à (^n-i) correspondant à W. 

Il s'ensuit que {A^-i) ^ (Bn) et on conclut par récurrence. 

4.4. Le lemme étant vérifié, l'assertion 4.2 (A) est 4.3 {A^) pour la matrice W égale à c{n + l)i-2 
avec la dernière colonne enlevée. 



§ 5. Formule (B) 

5.1. Maintenant on s'occupe de la formule 

P := {c{n)i + c(n) •'} • c(n) - c(n)ic(n)i_i = c(n - l)c(n + l)i_2 . 



On introduit n vecteurs de dimension n — 1, wi, . . . ,Wn qui sont les lignes de la matrice c(n + l)i-2 
sans la dernière colonne : 



/ è(l)2 
Kl)3 



VF 



Kl)3 

6(2)4 



6(2), 
6(1)„ b{2)n+l 
V 6(l)n+î-l 6(2)^+2 



6(l)n-l 
6(2)n 



6(l)n 
6(2)n+l 



6(ra - 2)2„-4 6(ra - 2)2„-3 
b{n - 2)2n-3 6(« - l)2n-2 

b{n - 2)2n+i-4 6(n - l)2n+i-3 y 



et n mineurs 



Wi = det{wi, . . . ,Wi, . . . ,w„ 



l,...,n . 



Par exemple, Wn = c(n), VF^-i = c(n)j_i, VF^-2 = c{n)'l. Donc, 

c(n + l)i_2 = b(n)2n+i-2Wn - b{n - l)2n-lWn-l 

+b{n - 2)2„_2W„_2 - . . . + (-l)"-i6(l)„+it^i 

= 6(n)2n+i-2W'„ - 6(n - l)2n-lWn-l + R 

OVi 

R = b{n- 2)2n-2Wn-2 - b{n - 3)2„-3W„_3 + . . . + {-ir-^b{l)n+lWi . 



(5.1.1) 



(5.1.2) 



5.2. On a n — 1 relations linéaires entre les Wi : la z-ième est obtenue en ajoutant à sa z-ième 
colonne et en développant le déterminant = par rapport à la dernière colonne. 
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Explicitement : 

b{n - l)2n+i-3Wn - b{n - l)2n-2Wn-l + 6(n - 2)2„-3t^n-2 - • • • + = , 

b{n - 2)2n+i-4Wn - b{n - 2)2n-3Wn-l + b{n - 2)2n-AWn-2 - • • • + {-XT'^biX^-lW^ = , 

b{2)n+iWn - b{2)n+lWn-l + b{2)nWn-2 - • • • + (-l)"-'6(2)4 ^2 + {-\T-''b{l)^Wi = , 
bil)n+i-lWn - bil)nWn-l + 5(l)„_iW^„_2 - . . • + {-ir-%{l)3W2 + (-l)"-^Kl)2^1 = . 



5.3. D'autre part, rappelons la matrice c(n)i : 



c(n)i = det 



/ Kl)2 6(1)3 

6(1)3 6(2)4 

6(l)„_i 6(2)„ 

V 6(l)„+i 6(2)„+2 



6(l)n-l 

b{2)n 



6(l)n \ 

6(2)„+i 



b{n - 2)2n-4 b{n - 2)2n-3 
b{n - 2)2n-2 b{n - l)2n-l j 



On développe cette quantité par rapport à la dernière colonne : 

c(n)i = b{n - l)2n-ic(n - 1) - 6(n - 2)2„-3M„_2 + . . . + (-l)"-i6(2)„+iM2 + (-l)"6(l)„Mi . 
Après la multiplication par — c(n)j_i = — Wn-\ on obtient : 

-c(n)ic(n)i_i = -5(n - l)2„-ic(n - \)Wn-\ (*) + -R' 

Oll 

i?' = b{n - 2)2„_3^n-lM„_2 - 6(n - 3)2n-4Wn-lM„_3 + . . . 

+(-l)"6(2)„+iH^„_iM2 + (-l)"-i6(l)„M/„_iMi . 



5.4. Maintenant remplaçons dans i?' les termes (— l)*6(n — ï)2n-i-^n-\ en utilisant les relations 



5.2 : 



6(n - 2)2„-3W^n-l = h{n - 2)2n+i-4Wn + b{n - 2)2n-4Wn-2 - ■ ■ ■ + (-l)'^-'6(l)n-lH^l , 



bi2)n+lWn-l = b{2)n+iWn + 6(2)„t^„_2 - • ■ ■ + (- 1)"-'6(2)4 1^2 + (-l)"-^6(l)3t^l , 
b{l)nWn-l = b{l)n+i-lWn + b{l)n-lWn-2 - • • • + (-l)"-26(l)3t^2 + (-1)^^-^6(1)2^^! 

Alors on obtient : 

-c{n)ic{n)i-i = -b{n - l)2n-ic(n - l)W„_i (*) 
+{6(n - 2)2„+i-4M„_2 - . . . + (-1)"6(2)„+,M2 + (-l)'^+^6(l)„+,_iMi} • c(n) + R", 

où : 

R" = |6(n - 2)2„-4W-„_2 - . . . + (-l)"-^6(l)„_iM/i| • M„_2 - . . . 
+(-1)" • |6(2),T^,_2 - . . . + i-ir-H{2)^W2 + (-l)"-^6(l)3l'Fi| • M2 
+(-1)^-1 • |6(l)„_iW-„_2 - . . . + (-1)^-26(1)3^2 + (-l)"-^6(l)2W^i} • M, 



1 • 
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5.5. Lemme. R" = c{n — 1)R. 

Démonstration. On introduit les vecteurs de dimension n — 2 : 

M = ((-l)"+iMi, (-1)"+2M2, . . . , M„_2) 

et 

6= (6(l)n+l,Kl)n+2,...,Kl)2n-2) • 

Alors la définition de R" se récrit : 

R" = M- C{n - 1) • W* (5.5.1) 

(où c(n — 1) = det C(n — 1), la matrice C(n — 2) étant symétrique) ; de plus, 

R = b-W* . 

Maintenant développons les quantités Mi par rapport à la dernière ligne : 

Mi = 6(l)2n-2Mi,„_2 - 6(l)2„-3Mi,„_3 + . . . + (-l)"+26(l)„+2Mi2 + (-l)'*+^6(l)„+iMa 

n-2 

= J2 i-ir^'bil)n+jMi„ i = 1, . . . ,n - 2 . 

On remarque que les quantités Mij sont les mineurs de la matrice (n — 2) x (n — 2) C(n — 1). Il vient : 

M = b-C{n-1) 

où 

c{n-i) = {{-iy+mij), 

donc C{n — 1) • C(n — 1) = c(n — 1). En substituant dans (5.5.1) : 

R" = b- C{n - 1) • C{n - 1) • W* = c(n - 1) • 6 • W* = c{n - 1)R, 

cqfd. 

5.6. Il s'ensuit que pour vérifier l'identité (B) il reste à démontrer que 

c{n)i + c{n)'l + b{n - 2)2„+i-4M„_2 - . . . 

(5.6.1) 

+(-l)'*6(2)„+iM2 + (-l)"+i6(l)„+i_iMi = b{n)2n+i-2c{n - 1) . 

Par contre, la quantité 

bin)2n+i-2cin - 1) - 6(n - 2)2„+i-4M„_2 + . . . + (-l)'^+^6(2)„+iM2 + (-l)'*6(l)„+i_iMi 
n'est autre que le développement de c(n)^ suivant la dernière colonne, donc (5.6.1) est équivalent à 

c{n)i + c(n)f = c{n% (5.6.2) 
qui a été déjà prouvée, cf. Corollaire 2.9. 

Ceci achève la démonstration du théorème 4.2, et donc du 1.5. 

5.7. Démonstration du théorème 1.9. En fait, nous l'avons déjà montré : la démonstration de la 
récurrence principale (3.3.2) n'utilise que les relations dans l'algèbre *B. 

Ces relations sont vérifiées si l'on définit les variables b{i)j à partir de coefficients de polynômes 
fi{x) et f2{x) comme dans 1.8, d'oti l'assertion. 
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DEUXIÈME PARTIE. 



POLYNOMES D'EULER ET DETERMINANT DE CAUCHY 

§ 1. Nombres 

1.1. Rappelons que pour un polynôme 

f{x) = a„x" + an-ix^'^ + . . . + ao 

les nombres h{j)i sont définis par 

i-i 

K3)i = (i- 2p)an-pan-i+p - j{n - i + j)an-jan+j-i • 

p=0 

On introduit les quantités : 



Ti : 



ai 

Qi-i aiai-2 



1i «1-1 



puis 



1.2. Par exemple : 



2n a„an_2 , 2n 

/3 1)2 = r • ^ 1 = • r„ - 1 , 

n - 1 <_i - 1 

o(-t\ (^nOin—i -, 

P(,J-jî 



On remarque que 

On définit les quantités 



n-i + 1 Qn-ian-i+l 
^n^n—i Qn—i+1 



«n-l«n-î+l Qn 



'^n—i+2fn—i+l ■ ■ - ^n ■ 



p=l 



(donc = 1 si i > j). Il s'ensuit : 



71 1 

n — i + 1 



1.3. De même : 

^n^n—i 0,n(^n—i 'ïn— l«n— i+l 

0^n-20.n-i-\-2 0'n-lO,n-i-\-l 0-n-20'n-i+2 

Par exemple : 



il}{n — i + 2, n)ip{n — i + 3,n 



"n-2 



tjj(n - 2, n)ip{n - 1, n - 1) = r„_2r^_ir„ 
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Il en découle 



^(2)- - ni anQn-i _^ n(i - 2) a„_iaTt-t+i _ ^ 

* n - i + 2 an-2an-i+2 n — i + 2 an-2an-i+2 

- i){n -i + 2, n)'4){n - z + 3, n - 1) H ^- — -^^{n - i + 3, n - 1) - 2 . 



n — i + 2 n — i + 2 

1.4. Un autre exemple : 

"''2''"'^ = V'(n - 4, n)V'(ra - 3, n - l)V'(ra - 2, n - 2) = r„_4r^_3r^_2r'^-irn • 



1.5. En général on pose : 

^{n,j,i) :-- 

et l'on aura : 



i-1 



- î + j + 9, « - 



i4 n(i - 2p) .... 
PU)î = y. ■(p{n-p,j -p,i-p) -J . 

p=0 

1.6. Passons maintenant aux déterminants c{n). On commence par un exemple : 



c(4) = det 



/ 6(1)2 6(1)3 6(1)4 
6(1)3 6(2)4 6(2)5 
\^ 6(1)4 6(2)5 6(3)6 



det 



^ {n - l)al_i(3{l)2 (n - 2)a„_ia„_2/3(l)3 (n - 3)a„_ia„_3/3(l)4 
(n - 2)an-ian-2/3(l)3 (^ - 2)a^_2/3(2)4 (n - 3)an-2an-3/9(2)5 
y (n - 3)a„_io„_3/?(l)4 (n - 3)an-2an-3/3(2)5 {n - 3)a^_3/?(3)6 

(n-l)/3(l)2 (n-2)/3(l)3 (n - 3)/3(l)4 
= (an-ian_2a„_3)' • det 1 (n - 2)/3(l)3 (n - 2)/3(2)4 (n - 3)/3(2)5 



(n-3)/3(l)4 (n-3)/?(2)5 (n - 3)/?(3)6 



1.7. En général 



c(m + 1) = [ n «n-j ] 



X det 



/ (n-l)/3(l)2 
(n-2)/3(l)3 



(n-2)/3(l)3 ... (n - m)/3(l)™+i \ 
(n-2)/3(2)4 ... (n - m)/3(2)^+2 



\ (n - m)/3(l)^+i {n-m)(3(2)m+2 



{n-m)/3{m)2m J 
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§ 2. Polynômes d'Euler et fonction hypergéométrique 



2.1. Suivant [Euler], on définit les polynômes 

1 



Enix) = -{(1 + ix/2n)^" + (1 - ix/2n)""} . (2.1.1) 

Donc, En{x) est un polynôme de degré 2n, avec le terme constant 1, ne contenant que des puissances 
paires de x. Plus précisément, 

£„(x) = |:(-l)'(2:)(|^- ("-2) 

Par exemple : 

Ei{x) = 1- , 

' 16 2048 65536 16777216 

2.2. Rappelons que la fonction hypergéométrique de Gauss est définie par 

^'^'^'^"^1-7 l-2-7(7+l) l-2-3-7(7 + l)(7 + 2) 

CX) 

= ^Cj(a,^,7)a;', 

1=0 

où 

a(a + l)...(a + z-l)-/3(/3+l)...(/3 + i-l) 
z!-7(7 + l)---(7 + ^-l) 

cf. [Gauss]. Il s'ensuit : 

Ci(-n/2,-n/2 + l/2,l/2) 

_ (-n/2)(-n/2 + 1) . . . (-n/2 + i - l) ■ (-n/2 + l/2)(-n/2 + 3/2) . . . (-ra/2 + i - 1/2) 
~ i!-(l/2)(l/2 + l)...(l/2 + i-l) 

_ (-l)*2-*n(n - 2) . . . (n - 2i + 2) • (-l)^2-*(n - l)(n - 3) . . . (ra - 2i + 1) 
~ i!-2-*-l-3-5...(2i-l) 

_ 2-^ • n(n - l)(n - 2) . . . (n - 2i + 1) _ f n\ 
~ 2-*-2-4...2M-3-5...(2i-l) ~ V2iy 

Donc 

F(-n/2, -n/2 + 1/2, 1/2,^^) = E [s^'i J = 2^^^ + ""^^ + " "^^"^ • ^^■^■^^ 

Il en découle : 

rF(-n/2, -n/2 + 1/2, 1/2, u'' /t") = ^{{t + u)" + (t - t^)"}, (2.2.2) 
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cf. [Gauss], no. 5, formula II. 

2.3. La formule (2.2.1) implique : 

= F{-n, -n + 1/2, 1/2, -x'^/An'^) 



(2.3.1) 



2.4. Si l'on écrit 

En{x) = Ç e„fe^2^ e^k := (-1)' ^fc) (2^ 



fc=o 

alors 



(2A;)!(2n)2'= (2fe)! V 2ny V 2ny ' ' ' V 2n 

d'où 

i.e. 

comme il faut. En d'autres termes, 

lim F(-n, -n + 1/2, 1/2, -x'^/An^) = cos x, 
OU, comme aurait pu écrire Gauss, 

F(-A;, A; + 1/2, 1/2, -xV4A;2) = cosx, 
k étant "un nombre infiniment grand" {dénotante k numerum infinité magnum). En fait, Gauss écrivit 

F{k, k', 1/2, -x^/4kk') =cosx, 
dénotante k,k' numéros infinité magnos, cf. [Gauss], no. 5, formula XII. 



§ 3. Asymptotiques 

3.1. On pose : 

/*)-Ê(-)'(^:)^^|:4"v. ,3.1.) 

Donc 

On désigne les quantités b{j)i,ri, etc. qui correspondent au polynôme en ajoutant l'indice (n) en 
haut : b{j)\"\r^^\ etc. 
Donc on aura : 
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X det 



(n - 2)mt^ in - 2)/3(2) 



3 
(n) 



(n) 



{n - m)/3(l);„ti \ 
(n-m)/3(2)^|2 

(n-m)/3(m)g / 



3.2. On a : 



= (-1) 



l 2î 



d'où 



,(n) _ af>a^'!>, _ [(2i-2)!]2[(2n-2i+2)!]^ 
î „{")2 (2i) ! (2n-2i)! (21-4)! (2n-2î+4)! 

1 

_ (2i-2)(2i-3) (2ra-2i+l)(2n-2i+2) 
~ 2i(2V-l) ■ (2n-2i+3)(2n-2i+4) ■ 



En remplaçant i par n — i, 



(„) _ (2î + l)(2i + 2) (2n - 2i - 2)(2n - 2i - 3) 
^"-^ ~ (2i + 3)(2i + 4) ' (2n-2i)(2ra-2i-l) 

On s'intéresse aux valeurs limites : 



(n) ^ (2i + l)(2i + 2) 
'^-^ (2z + 3)(2i + 4) ■ 



Il s'ensuit : 



ip{(Xi — i + 2, oo) := lim ip(n — i + 2,n) = 



1 • 2 



(2z - l)2i 



^{oo — i + 3, oo — 1) 
■^(oo — i + 4, oo — 2) 



3 ■ 1 



(2i-2)(2î-3) 

5 • 6 
(2z-4)(2i-5) 



etc. 



3.3. Maintenant on veut calculer 

Il est commode de poser : 
On a : 

d'où 



B{j)r ■■= mt^ + j ■ 



=i-V(oo-i + 2,oo) 



2z- 1 



(oo) _ 2(i - 1) 

2i-\ 



Ensuite, 



B{2)'f°^ = i ■ il){oo -i + 2, oo)'(/'(oo - i + 3, oo - 1) + (i - 2) • ^(oo - i + 3, oo - 1) 



= V(oo-z + 3,oo-l)-|s(l);°°)+i-2| = 



3-4 



(2i-2)(2i-3) [2i-l 



1 



+ i-2 



3 • 2 

2i-l' 
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d'où 



m 



(oo) _ 4(i - 2) 



2i - 1 



De même, 
5(3)S°°^ = ij{oo - z + 4, oo - 2) ■ |s(2) 



+ Z-4 



5 • 6 



d'où 



(oo) 



3 • 2 

(2i-4)(2i-5) l2i - 1 

6(i - 3) 
2z - 1 ■ 



+ i - 4 



3.4. En général, la récurrence évidente fournit 

s(oo) _ (2j - 1) • j 



et 



2i-l 

(oo) ^ 2j(/ - j) 
^ 2i-l 



3.5. On définit les nombres 



c(m + 1)°° := det 



/m+l 



V/3(l)ïïi /?(2) 



/3(2)f) 



(CX)) 

m+2 



/?(2)ÏÏ2 



Donc on aura : 



n 



i=l 



• c(m + = c(m + 1)~ • n"" + 0(n"*-^) . 



Les calculs précédents fournissent par exemple : 



c(4)°° = det 



2 _4 _6 

"3 5 7 

4 _8 _ 12 

■5 7 9 

6 12 18 



(_1)3 . 2 . 4 • 6 • det 



\ D _io / 

\ 7 9 11 / 

= (-1)3-23 . (3!)2 - det 



/ 1 2 3 

/ 3 5 7 

12 3 

5 7 9 

\ i 2 J. / 

\ 7 9 11 / 



/ 1 1 1 
/ 3 5 7 

111 
5 7 9 



1 1 X y 
9 11 ' 



En général on obtient 



c(m + 1)°° = (-1)™ • 2™ • (m!)2 • det 



1 
5 



2m+l 
1 

2m+3 



\ 2m+l 2m+3 



5 • 3 

2i- 1^ 



• 4m-l / 

On remarque que la dernière matrice (une variante de la matrice de Hilbert) est du type Hankel. 
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3.6. Le déterminant 



C{m + 1) := det 



1 

5 



\ 2m+l 2m+3 



2m+ 
1 

2m+3 



4m- 1 / 



est un cas particulier du déterminant calculé par Cauchy (d'oii le caractère £), cf. son Mémoire sur les 
fonctions alternées et sur les sommes alternées, pp. 173 - 182 dans [Cauchy]. 

Rappelons que, étant données deux suites xi, . . . , et yi, ■ ■ ■ ,ymi le théorème de Cauchy dit que 



det((xi + %) ^)T,j=i 



d'où, en posant Xi = 2i — 2, j/j = 2i + 1, 

<L{m + 1) 



^ Y\l<i<3<m jxj - Xi){yj - Vi) 

ni}=i (^î + ^j) 
n^,=i"(2i + 2j-i) • 
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